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ABSTRACT 

On the basis of a recently derived color-metallicity relation and stellar parameters from the Sloan 
Digtal Sky Survey Data Release 7 spectroscopic survey, a large sample of red horizontal branch 
candidates have been selected to serve as standard candles. The metallicity and kinematic distributions 
of these stars indicate that they mainly originate from the thick-disk and the halo populations. The 
typical thick disk is characterized by the first group peaking at [Fe/H] ^ —0.6, V-^ot ^ 170kms~^ 
with a vertical scale height around \Z\ ^ 1.2 kpc, while stars with [Fe/H] < —0.9 are dominated by 
the halo population. Two sub-populations of the halo are suggested by the RHB stars peaking at 
[Fe/H] ^ —1.3: one component with V^ot > Okms~^ (Halo I) shows a sign of metallicity gradient in 
the [Fe/H] versus \Z\ diagram, while the other with V^ot < Okms"^ (Halo II) does not. The Halo I 
mainly clumps at the inner halo with i? < 10 kpc and the Halo II comes both from the inner halo 
with i? < 10 kpc and the outer halo with R> IQ kpc based on the star distribution in the R versus 
\Z\ diagram. 

Subject headings: stars:horizontal-branch 
the halo 



Galaxy: kinematics - Galaxy: the thick disk - Galaxy: 



1. INTRODUCTION 

Red horizontal-branch (RHB) stars are found in 
many Galactic globular clusters (GCs) with intermedi- 
ate metallicity —1.7 < [Fe/H] < —0.3 (see, e.g., Piotto 
et al. 2002). This metallicity range includes stars from 
the halo, the thick disk, the thin disk, and probably the 
bulge. Due to their nearly constant absolute magnitude, 
this type of stars constitutes an ideal sample to trace 
various Galactic populations in the Galaxy. 

Generally, the identification of RHB stars in the field 
is rather difficult just by their colors but some attempts 
to identify this type of star had been made (Rose 1985; 
Tautvaisiene 1996; Preston 2006). However, these works 
did not produce a large sample of RHB stars (for ex- 
ample, a sample with star number larger than 100) for 
statistic work. With the aid of Hipparcos parallaxes, 
Kaempf et al. (2005) selected a large sample of RHB 
stars in the field and calculated their kinematical pa- 
rameters based on available radial velocities and proper 
motions. They found two populations: a disk one having 
\Z\ ~ 0.6 kpc and a halo one showing \Z\ ~ 4 kpc. This 
is the only work to trace the structure of the Galaxy via 
a large sample of RHB stars. Due to the distance limi- 
tation of the Hipparcos survey, these results were based 
on stars in the solar neighborhood, where the majority 
of the disk stars are from the thin disk with a minority 
from the thick disk and only a sprinkling of halo stars. It 
is of high interest to trace the disk and halo populations 
via RHB stars locating far away from the solar neighbor- 
hood and well above the Galactic plane. For this pur- 
pose, planned or ongoing kinematically unbiased surveys. 
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such as the Global Astrometric Interferometer for Astro- 
physics (GAIA; Turon et al. 2005), the SDSS/SEGUE 
(York et al. 2000; Yanny et al. 2009), and the LAMOST 
(Zhao et al. 2006) projects, will provide a chance to ob- 
tain and to study a large sample of RHB stars to the 
distant halo. 

In this work, we attempt to investigate the properties 
of RHB stars and to carry out kinematic traces of various 
Galactic populations using RHB stars selected from the 
SDSS spectroscopic survey. As compared with the previ- 
ous study by Kaempf et al. (2005), the present work has 
several advantages. First of all, the SDSS survey covers 
quite distant regions of the Galaxy, and our study can be 
extended far away from the solar neighborhood. Second, 
spectroscopic data of the SDSS survey provide stellar pa- 
rameters, radial velocities, and metallicities. This data 
set not only makes it possible to identify RHB candi- 
dates more easily, but also provides the most powerful 
tool to investigate the Galactic populations by a combi- 
nation of chemical and kinematical ingredients, in view 
of the fact that our Galaxy evolves both chemically and 
dynamically. Finally, it is expected that RHB stars se- 
lected from the SDSS survey could extend to a much 
lower metallicity range than that of Kaempf et al. (2005), 
which mainly concentrated on solar metallicity, and thus 
they will contribute much to the understanding of the 
thick-disk and the halo populations. 

2. SELECTION OF RHB STARS 

RHB stars are easily seen based on the color-magnitude 
diagrams (CMDs) of GCs, and statistically it is possible 
to pick out this type of star in the field of the solar neigh- 
borhood with the aid of Hipparcos parallaxes. Without 
previously-known distance information for stars far from 
the solar neighborhood, the selection of RHB stars in the 
field is rather difficult. But the problem of identification 
is much easier if the metallicity of the star is known, as 
shown in Straizys et al. (1981). Based on high resolution 
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Figure 1. {g — r)o-[Fe/H] diagram for clusters in Chen et al. 
(2009). The solid line indicates a linear fitting to the data. 

and high signal-to-noise ratio spectra, around 30 RHB 
stars have been identified in the literature (e.g., Behr 
et al. 2003; Carney et al. 2003, 2008). In order to en- 
large this sample, RHB stars have been identified from 
the SDSS low resolution (R ~ 1800) spectroscopic sur- 
vey. The SDSS Data Release 7 (DR7; Abazajian et al. 
2009) provides photometry and spectroscopic mctallici- 
ties for a large sample of high-latitude stars, allowing the 
identification of a large sample of RHB stars in the field. 

2.1. The Metallicity- Color Relation 

As noticed in Chen et al. (2009), the {g — r)o values of 
RHB stars in metal-poor GCs with [Fe/H] = —1.6 dex, 
such as Pal 3, NGC 7006, and M3, are about 0.25-0.30 

mag and extend to the red, reaching 0.60 in M71 with 
[Fe/H] = -0.8 and to 0.98 in NGC 6791 with [Fc/H] = 
-1-0.4. In fact, there is a correlation between {g — r)o and 
[Fe/H] for RHB stars in clusters. In order to obtain a 
quantitative relation, we perform a linear regression to 
the data in Chen et al. (2009), which gives 
(g - r)o = 0.343(±0.039)[Fe/H] -f- 0.829. 
The scatter around this relation is approximately 0.07 
mag. Figure 1 shows the peaking [g — r)o (and its error) 
versus [Fc/H] for each cluster and a linear relation is 
obtained by fitting to these points. 

Assuming that RHB stars in the field follow the same 
(gf — r)o and [Fe/H] relation as these clusters, the first 
step is to apply this color-mctallicity relation to all stars 
in SDSS-DR7 where spectroscopic metallicities and pho- 
tometric colors are available in the published catalog. 
In the present work, the interstellar reddening values 
E{B — V) are based on the extinction map of Schlegel 
et al. (1998) and we adopt the reddening corrections for 
ugriz colors from Fukugita et al. (1996). In the selection, 
a deviation in {g — r)o to the above relation 5 < 0.15 mag 
is used to select stars; this value is adopted by taking into 
account an error of 0.2 dex in [Fe/H] in SDSS-DR7, an 
error of 0.1-0.2 mag in the reddening uncertainty, and 



an approximate (g — r)o color extension of 0.1 mag for 
RHB stars in cluster. Then, we select stars with signal- 
to-noise ratios of their spectra larger than 10 and metal- 
licity larger than [Fe/H] > —2.0 so that their radial ve- 
locities and stellar parameters can be reliably determined 
from the SDSS-DR7 catalog. Finally, an upper cutoff 
at 5 = 20 mag is necessary because the proper motions 
from the USNO survey become seriously contaminated 
by misidentifications toward faint magnitudes according 
to Fuchs et al. (2009). But the effective cutoff in mag- 
nitude of our final sample is around 5 = 18 — 18.5 mag 
and the saturation flag in the SDSS-DR7 catalog is used 
to avoid stars whose magnitudes suffer from saturation 
effects. With these criteria, we obtain our preliminary 
sample of ~ 120, 000 stars. 

2.2. The Toff versus logg diagram. 

Using the above {g — r)o [Fe/H] relation, blue hori- 
zontal branch (BHB) and RR Lyrae stars are excluded 
at the bluer colors, and red giant branch stars at the 
red limit, for a given metallicity. However, the sample is 
significantly contaminated by main-sequence stars (MS), 
turnoff (TO), and subgiant (SG) stars with the same 
metallicity; these can be distinguished if their luminos- 
ity or gravities are known. In the CMD of the clus- 
ters in Chen et al. (2009), at a given metallicity, RHB 
stars have brighter absolute magnitudes than those of 
AIS/TO/SGB stars despite their similar colors, and the 
separation between RHB and MS/TO/SGB stars in ab- 
solute magnitude at a given color becomes larger as the 
metallicity decreases. In view of this, stellar parameters, 
i.e., temperature, gravity and metallicity, taken from the 
SDSS-DR7 catalog, can be used to identify RHB stars. 
The Tcff— logg diagram for the first 20000 stars in the pre- 
liminary sample and the corresponding contour map is 
shown in Figure 2. This figure shows that MS/TO/SG 
and RHB stars can be distinguished with RHB stars 
clumping over the range of 4500-5900 K in temperature 
and 1.8 3.5 dex in logg, which are the selection criteria 
for our sample stars. In order to check these tempera- 
ture and gravity ranges for RHB stars, we have compiled 
stellar parameters of previously identified RHB stars in 
the field and in clusters based on high resolution spec- 
tral analyses in the literature. Most works (Behr et al. 
2003; Carney et al. 2003, 2008) agree with the ranges of 
4600 6200K in temperature and of 2.0 3.3 dex in ?ogg for 
RHB stars, which are consistent with the defined ranges 
of RHB stars in Figure 2. 

3. STELLAR DISTANCES AND SPACE 

VELOCITIES 

3.1. Distances, proper motions, and radial velocity 

Stellar distances have been calculated rising the 
absolute magnitude-metallicity relation of Afg = 
0.492[Fe/H] + 1.39 derived by Chen et al. (2009). We 
have compared them with distances derived from the i 
magnitude and Mf — [Fe/H] calibration derived in Chen et 
al. (2009), and they give similar distances. The radial ve- 
locities and proper motions are based on values provided 
by the SDSS-DR7 catalog. Note that the proper motions 
in the SDSS-DR7 catalog have been corrected for the sys- 
tematic error noticed by Munn et al. (2008). Finally, only 
stars with high quality data in photometry {ug < 0.1, 
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Figure 2. Tgg - logg diagram and its contour map for the first 
20000 stars in the preliminary sample. The box indicates our se- 
lected RHB candidates. 

ai < 0.1 rnag), in proper motions {a^ < 3mas yr~^), 
and in radial velocity {(Jrv < lOfcm s~^) are included in 
the sample. 

3.2. Spatial velocity 

Stellar positions (X,Y,Z) and heliocentric space veloc- 
ity (U,V,W) are calculated following the same procedures 
as in Chen ct al. (2000). Wc use a left-handed system 
in which U is positive toward the Galactic anti-center. 
The solar motion of (U,V,W)0 = (7.5, 13.5, 6.8) km g-^ 
with respect to the LSR from Francis & Anderson (2009) 
and the local standard of rest (LSR) velocity of 220 km/s 
(corresponding to Vhei ^ —225km s~^) are adopted in 
the present work. 

The errors in spatial velocity come from the uncertain- 
ties of radial velocity, proper motion and distance. The 
error in radial velocity of 3 km/s corresponds to a spa- 
tial velocity error of 5 km, s^^. The distance error of 
10%, estimated from an error in absolute magnitude of 
0.2 mag, corresponds to a velocity error of 5-15 km, s~^. 
The largest error comes from uncertain proper motions. 
As shown in Ivezic et al. (2008), a random error of 3 
m,as yr~^ in proper motions will lead to a velocity er- 
ror of about 15 kms~^ at a distance of 1 kpc and about 
80 km at a distance of 5 kpc. Although the total 
error in spatial velocity is quite large, the data still rep- 
resent usable measurements for a large sample because 



the systematic errors are much smaller (^10 km at 
a distance of 7 kpc) according to Ivezic et al. (2008). 

In principle, a total velocity larger than 600kms~^ is 
unreasonable, or indicates stars unbound to our Galaxy. 
Such stars are usually excluded from the final sample 
(e.g., Klcmcnt ct al. 2009). In the present work, we 
keep them in our sample for further check, but we have 
checked that the main results arc generally the same 
when we excluded them from the sample. With this 
in mind, we avoid to overexplain the results based on 
the spatial velocities only and instead statistical results 
based on mctallicity and distance to the Galactic plane 
{\Z\) are more favorable in the following analysis. 

The final sample includes 5391 stars for further anal- 
yses. Figure 3 shows the distributions of distances, g 
magnitude, reddening, and galactic latitudes of the fi- 
nal sample. Their identifications, stellar parameters, and 
spatial velocities are published electronically and a sam- 
ple table consisting of the first 10 stars is presented in 
Table 1. 

4. RESULTS AND DISCUSSIONS 

4.1. The metallicity and kinematical distributions 

By choosing a single stellar type, RHB stars, it is in- 
teresting to investigate the metallicity distribution and 
compare our results with those from other types of stars, 
e.g., MS or TO stars. In Figure 4, a histogram of our 
sample shows two peaks in the metallicity distribution 
with a division at [Fe/H] ~ —0.9 dex. Two Gaussian fits 
to the distribution show a metal-poor sub-population at 
Fe/H] ~ —1.3 dex, and a metal-mild sub-population at 
Fe/H] ~ —0.6 dex. As compared with the result from 
MS stars within 4 kpc of the Sun in the SDSS I and II sur- 
vey by Carollo et al. (2010), we find that the two peaks in 
the metallicity distribution are quite similar. They sug- 
gested that stars in the metal-mild peak belong to the 
disk population (mainly the thick disk), while stars in 
the metal-poor peak belong to the halo population. How- 
ever, there are some diff'erences in star numbers for each 
population between RHB stars in the present work and 
MS stars in Carollo et al. (2010). In our work, the metal- 
poor component, peaking at [Fc/H] ^ —1.3, outnumbers 
the metal-mild component, peaking at [Fe/H] ~ —0.6, 
which is opposite to the result in Carollo et al. (2010), 
where the metal-mild component surpasses the metal- 
poor component. This is easily understood considering 
that their sample of stars is limited to d < 4 kpc from 
the Sun, and thus the disk population contributes signif- 
icantly to their sample, while RHB stars in our sample 
have brighter absolute magnitudes and extend to more 
distant regions of the Galaxy where the halo dominates. 

The variations in the metallicity distribution with dif- 
ferent distances to the Galactic plane {\Z\) are shown 
in Figure 5, where the metallicity peak shifts from 

[Fe/H] 0.4 dex for \Z\ < 1.5, '[Fe/H] 0.6 dex 

for 1.5 <\Z\ <5 to [Fe/H] 1.3 dex for \Z\ > 8. At 

5 < \Z\ < 8, the star numbers between the two compo- 
nents are comparable. Note that the result in Figure 20 
of Carollo et al. (2010) shows the shift of the metallicity 
peak from —0.6 at kpc < \Z\ < 1 kpc, to —1.3 at 3 kpc 
< jZj < 4 kpc, reaching to —1.6 at 6 kpc< \Z\ < 7 kpc 
and to -2.2 at \Z\ > 9 kpc. Therefore, the exact [Fe/H] 
peaks at a given \Z\ between Carollo et al. (2010) and 
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Table 1 

Identifications, stellcir parameters, coordinates, and spatial velocities of the first ten stars in the final sample are presented. 
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the present work are somewhat different due to the selec- 
tion of different spectral types of stars in the sense that 
MS/TO stars in their sample include significant contribu- 
tions from the thin disk and the halo with [Fe/H] < —1.6, 
while our sample of RHB stars has the largest contribu- 
tion from the thick disk, the minority contribution from 
the halo and nearly negligible contribution from the thin 
disk. 

In order to investigate the origins of the two popu- 
lations in the metallicity distribution for RHB stars in 
the present work, the Toomre diagrams for stars with 
[Fe/H] > —0.9 and [Fc/H] < —0.9 arc shown in Figure 6 
where the solid circle indicates the division between the 
thick disk and the halo with Vtot = 210 km s"""^ following 
Nissen & Schuster (2010). Note that there is a signif- 
icant number of stars with a total velocity larger than 
600kms~^ in these figures due to the very large errors in 



the determination of space velocities when the distance 
of the star is larger than 10 kpc. Therefore, we do not 
attempt to draw any firm conclusion from the kinemat- 
ical data, but the statistical argument for dividing stars 
into different groups may be right. 

It is clear that the metal- mild component of RHB stars 
clumps within the thick-disk region with some extension 
to the thin-disk region, while the majority of stars in 
the metal-poor component are located in the halo region. 
This result is consistent with that of Carollo ct al. (2010) 
in the sense that the metal- mild component, peaking at 
[Fc/H] ^ —0.6, mainly originates from the thick disk and 
the metal-poor component, peaking at [Fe/H] ~ —1.3, 
generally belongs to the halo. In Figure 7, it is inter- 
esting to compare the kinematic diagrams between the 
two metallicity peaks for stars with 5 kpc < |Z| < 8 
kpc in Figure 6 (the third panel). Again, it shows that 
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Figure 4. Mctallicity distributions of RHB stars. The dotted 
lines show the Gaussian fits to the two populations separated by 
the metallicity division at [Fe/H] ~ —0.9. 
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Figure 5. Metallicity distributions of RHB stars at different \Z\. 



the metal-poor component at [Fe/H] ~ —1.3 mainly has 
the halo kinematics and the metal-mild component at 
[Fe/H] ~ —0.6 has the thick disk kinematics. In the fol- 
lowing sections, wc separate the stars into two parts with 
[Fe/H] > -0.9 and [Fe/H] < -0.9, respectively 

4.2. The [Fe/H] versus \Z\ diagram 

Figurere 8 shows the distribution of 5391 RHB stars 
in the [Fe/H] versus \Z\ diagram with the pink color 
corresponding to a high density, green to middle den- 
sity, and black to the lowest density. The division be- 
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Figure 6. Toomrc diagrams of RHB stars for [Fe/H] > —0.9 and 
[Fe/H] < —0.9. The circles indicate the division between the thick 
disk and the halo following Nissen & Schuster (2010). 



tween the thick disk and the halo seems to be quite clear 
in the [Fe/H] versus \Z\ diagram as clearly shown in 
Figure 8, where the edge of the thick disk could be as 
high as \Z\ ~ 8 kpc at [Fe/H] ~ —0.6 and reduces to 
|Z| ~ 2 kpc at [Fe/H] ^ —1.5. Figurere 9 (upper panel) 
shows individual stars in this diagram. It seems that 
there is a metallicity gradient in the vertical direction of 
the Galaxy, but the scatter around the gradient is quite 
large. The \Z\ range increases with decreasing metallicity 
probably due to the overlapping metallicity of stars from 
the halo, the thick-disk and the thin-disk populations at 
this metallicity range. Actually, the metallicity gradi- 
ent is more significant in the middle metallicity range of 
— 1.2 < [Fe/H] < —0.4 than the two outer ranges in this 
[Fe/H] versus \Z\ diagram. 

In order to investigate the properties of the thick- 
disk population, we try to exclude the contribution 
from th(^ halo population and plot in the lower panel 
of Figure 9 the [Fe/H] versus \Z\ diagram for stars with 
\Z\ < 5 kpc and [Fe/H] > —0.9, which mainly origi- 
nate from the thick disk. This figure shows that the 
trend in the [Fe/H] versus \Z\ diagram is quite sig- 
nificant and the scatter around the gradient is greatly 
reduced. A linear regression to the data indicates a 
relation of [Fe/H] = -0.292(±0.012)|Z| + 0.116 with 
a scatter of 0.037 dex. Taking into accomit the fact 
that a few points with [Fe/H] > —0.35 do not fol- 
low the above trend and may belong to the thin disk, 
they are excluded from the fitting, and the result is 
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Figure 7. Toomrc diagrams of RHB stars with 5 < \Z\ < 8 kpc 
for [Fe/H] > -0.9 and [Fe/H] < -0.9, respectively. The symbols 
are the same as Figure 6. 

[Fe/H] = -0.255(±0.014)|Z| + 0.020 with a scatter of 
0.023 dex. The six points with [Fe/H] > -0.35 follow 
the trend of [Fe/H] = -0.467(±0.083)|Z| + 0.323 with a 
scatter of 0.031 dex for \Z\ < 1.5 kpc. 

4.3. The V[ot versus \Z\ diagrams 
Figurere 10 shows the kinematical gradients, K-ot {— 
y + 220 kms~^) versus \Z\, and Figures 11 and 12 show 
the U and W velocities versus \Z\ for RHB stars; for 
both figures, [Fe/H] < -0.9 and [Fe/H] > -0.9 in 
the upper and lower panels, respectively. These fig- 
ures show that ([/) and (W) velocities generally do not 
vary much with \Z\, independent of metallicity, while for 
[Fe/H] < —0.9 the (K-ot) tends to decrease with increas- 
ing \Z\ for the region \Z\ < 10 kpc and becomes flat for 
the region \Z\ > 10 kpc. The rotational- velocity gradi- 
ents are approximately —30 km s^^ kpc^^ for stars with 
[Fe/H] < -0.9 and \Z\ < 10 kpc, and about -19 km a'^ 
kpc-i for stars with [Fe/H] > -0.9 and \Z\ < 15 kpc. 

Quantitatively, different authors give somewhat dif- 
ferent gradients for \Z\ < 5 kpc. For example, Car- 
ollo et al. (2010, in their Figure 11) found a decreasing 
(K*) with \Z\ for stars with -0.8 < [Fe/H] < -0.6 and 
\Z\ < A kpc with a rotational-lag gradient of —36 km 
s~^ kpc~^. Chiba & Beers (2000) found a rotational-lag 
gradient of —30 km s~^ kpc~^ based on a large sample 
of stars (- 1200) with \Z\ < 2 kpc. Allende Prieto et al. 
(2006) adopted SDSS photometric distances and radial 
velocities (without using proper motions), and found a 
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Figure 8. Distribution of 5391 RHB stars in the [Fe/H] vs. \Z\ 
diagram. The star density decreases as the color varies from the 
pink, the blue, the green to the black. 




Figure 9. [Fe/H] vs. \Z\ diagram for the full sample (upper panel) 
and for stars with [Fe/H] > —0.9 and \Z\ < 5 kpc (lower panel). 
The moan \Z\ at a given [Fe/H] and the scatter are shown by the 
large filled circles and solid lines, respectively 



rotational-lag gradient of —16 km s ^ kpc ^ for stars be- 
tween \Z\ — 1 and 3 kpc. In JVIajewski (1992), the value 
is —21 km s~^ kpc~^ based on deep proper-motion sur- 
vey of the Galaxy out to \Z\ 6 kpc. Our value of —19 
km s~^ kpc~^ for stars with [Fe/H] > -0.9 and \Z\ < 15 
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Figure 11. U velocities vs. \Z\ for stars with [Fe/H] < —0.9 and 
[Fe/H] > —0.9 (upper and lower panels, respectively). 



Figure 10. Vrot versus \Z\ for stars with [Fe/H] < 
[Fe/H] > —0.9 (upper and lower panels, respectively). 



-0.9 and 



kpc is quite consistent with previous values within the 
error range. 

It has been suggested that the rotational lag and the 
velocity dispersions vary with distance from the Galactic 
plane (Majewski 1994). However, our work shows that, 
for [Fe/H] < —0.9, the (Kot) decreases with increasing 
\Z\ for the region of \Z\ < 10 kpc, while for \Z\ > 10 
kpc, where the halo dominates, there is no significant 
trend. In Bond et al. (2010) the halo population with 
[Fe/H] < —1.1 does not show any trend in the (K-ot) 
with \Z\ over the region of jZj < 5 kpc, which does not 
agree with our result for [Fe/H] < —0.9 and \Z\ < 10 
kpc. This discrepancy may be easily understood if we 
assume that stars with [Fc/H] < —0.9 in the present 
work consist of two halo sub-populations. Finally, the 
scatters around these trends in Figure. 10-12 arc rather 
large, and further studies with larger samples and high- 
quality data are needed. 

4.4. The [Fe/H] vs. Vrot diagram 

As described above, this diagram is a useful way to 
trace the structure of the Galaxy when abundance and 
kinematical data are combined. Figure 13 shows the 
histograms of Vrot (upper panel) and the [Fc/H] versus 
V^ot diagram (lower panel) for RHB stars. The upper 
panel shows that the metal-mild component peaks at 
Vrot ~ 170kms~^ and the metal-poor component spans 
a wide range in V^-ot without any sharp peak. That is, 



the Vrot dispersion in the former is smaller than the lat- 
ter component. This is expected in the context of the 
Galactic stellar populations with the halo population 
with [Fe/H] < —1.0 having a large velocity dispersion 
as compared to the thick-disk population, with a mctal- 
licity peak at [Fe/H] ^ —0.6, and the thin-disk popula- 
tion with a mctallicity peak at [Fe/H] ^ —0.2. Actually, 
the distribution in Vrot for the metal-poor component is 
rather broad and in this work we separate the metal- 
poor component into two groups with Vrot > Okms""'^ 
and Vrot < Okms^^ in the lower panel of Figure 13. 
The metal-mild population with [Fc/H] ^ —0.6 has a 
peak at V^ot ~ 170kms~^, which is quite similar to that 
derived from the solar neighborhood by Soubiran et al. 
(2003), who derive a rotational lag (Strombcrg asymmet- 
rical drift) of Viag ^ 51kms^^ with respect to the LSR, 
and is somewhat lower than the average rotational lag 
of Viag ~ 80kms~^ by Fuhrmann (1998) for a group of 
16 thick-disk stars. For the metal-poor population, we 
separate them into two sub-populations with the divi- 
sion of stars with Vrot > Okms^^ being "Halo I" and 
stars with V-ot < Okms^^ being "Halo 11". In Figure 14, 
the [Fe/H] versus |Z| diagrams for the Halo I and II sub- 
populations are shown in the upper and lower panels, 
respectively. It shows that there is a hint of mctallicity 
gradient in the Halo I, as for the thick-disk population, 
while for the Halo II population the mctallicity gradient 
is much weaker, if it exists at all. This gives concordance 
to our separation of the halo population into the Halo I 
and the Halo II components for stars with [Fe/H] < —0.9 
in Figure 13. 
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Figure 12. W velocities vs. \Z\ for stars with [Fe/H] < —0.9 and 
[Fe/H] > —0.9 (upper and lower panels, respectively). 

Actually, the division of the halo at Viot ^ Okms"-^ is 
consistent with the result of Carollo et al. (2007), who 
found an inner halo with a slightly prograde rotation and 
a outer halo with a net retrograde. In their Figure 5, the 
division between the inner and the outer halo is quite 
clear in the R — \Z\ diagram. In view of this, it is inter- 
esting to investigate the division for the halo population 
based on RHB stars. Figure 15 shows the distribution of 
RHB stars for the Halo I and the Halo II sub-populations 
in the R—\Z\ diagram. The main result from this figure 
is that the Halo I stars clump within i? < 10 kpc and 
\Z\ < 10 kpc, while the Halo II stars show two clumps 
with one within 10 kpc and the other grouping outside 
10 kpc. In this sense, we favor the suggestion that the 
transition between the inner halo and the outer halo is 
around 10 kpc. Note that the Figure 5 of Carollo et al. 
(2007) is divided by the mctallicity while our Figure 15 
is separated by the rotation, and thus the comparison 
bc;twec>ii the two works is not direct. However, the struc- 
ture of the halo consisting the inner and the outer parts 
is the same. 

4.5. Implications for the evolution of the Galaxy 

Since the first suggestion for the thick disk in our 
Galaxy was made on the basis of star counts by Gilmore 
& Reid (1983), a considerable number of studies have 
helped to characterize the Galaxy's thick-disk compo- 
nent, in an effort to deduce its origin and better under- 
stand its nature. Most works suggest that the thick disk 
consists of old (> 10 Gyr) stars of intermediate metallic- 
ities, (-1.0 < [Fe/H] < -0.2) with a rotational lag of w 50 



Figure 13. Histograms of Vlsr (upper panel) and the [Fe/H] 
versus Vlsr diagram (lower panel). The thick-disk population 
consists of stars with [Fe/H] > —0.9, while the halo population 
is divided into two sub-populations, the Halo I (crosses) and the 
Halo H (dots). 

km s~^, and a scale height between 0.8 and 1.2 kpc. 
Our results from RHB stars with [Fe/H] < —0.9 dex are 

generally consistent with these numbers; this population 
peaks at [Fe/H] - -0.6±0.2 dex, V.^t ~ 170±20kms-\ 
and a vertical scale height of 1.2 ± 0.3 kpc. 

However, the thick disk may have a metal- weak com- 
ponent (MWTD) (e.g.. Beers et al. 2002). Based on 
MS/TO stars from the SDSS survey, Ivezic et al. (2008) 
reveal a feature consisting of stars of [Fc/H] —1.0 
dex and \Z\ ^ 2 — 3 kpc, which they classified as the 
MWTD population. Brown et al. (2008) presented a 
complete spectroscopic survey of 2414 Two Micron All 
Sky Survey (2MASS)-selected BHB candidates selected 
over 4300 dcg^ of the sky and found that the BHB stars 
located at a distance from the Galactic plane \Z\ < 4 kpc 
trace what is clearly a metal- weak thick-disk population, 
with a mean mctallicity of [Fe/H] = — 1.7 dex. But the 
existence of MWTD among BHB stars with \Z\ < 4 kpc 
is not confirmed by their further work in the paper by 
Kinman et al. (2009). 

In the present work, the MWTD has not been found as 
shown in Figure 14 (upper panel), where there are only 
a few stars with [Fe/H] < —0.9 and \Z\ < 5 kpc, and 
very few with \Z\ < 2 kpc. The Halo I sub-population 
with [Fe/H] ^ —1.3 and Vrot > Okms"-*^, generally has 
\Z\ > 5 kpc. Since it is well accepted that the edge of 
the thick disk is at about 5.5 kpc above the Galactic 
plane as suggested by Majewski (1994), and since we 
have measured a scale height for the thick disk of \Z\ = 
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Figure 14. [Fe/H] vs. \Z\ for the Halo I (upper panel) and II 
(lower panel) sub-populations. 

1.2 ± 0.3 kpc, we suggest this sub-population, Halo I, 
really corresponds to the halo. 

In connection with the halo population, our result con- 
firms with the existence of the inner and the outer halo 
sub-population in the R - \Z\ diagram. From the [Fe/H] 
- \Z\ diagram, it hints that the inner halo shows a metal- 
licity gradient while the outer halo does not. From our 
data, we cannot detect any sign of stellar streams re- 
ported before due to the selection effect of the spectro- 
scopic survey of the SDSS project and the star number 
in our sample is too small for detecting any overdensity 
produced by stellar streams. 

5. CONCLUSIONS 

Based on photometric and spectral data, a group of 
RHB stars have been selected from the SDSS survey, 
whose distances can be estimated by using the calibra- 
tions from Chen et al. (2009). Combining with available 
proper motions from the USNO survey and radial veloc- 
ities from the SDSS survey, the space velocities (U,V,W) 
have been calculated for a sample of 5391 RHB stars, and 
their metallicities and kinematical gradients have been 
investigated in order to trace the evolution of the Galaxy. 
The main results are as follows. (1) There are two peaks 
in the metallicity distributions of RHB stars with a di- 
vision at [Fe/H] 0.9. Stars with [Fe/H] > -0.9, 

centered at [Fe/H] ^ —0.6, correspond to the thick- 
disk population, while stars with [Fe/H] < —0.9, peak- 
ing at [Fe/H] ^ —1.3, have generally halo kinematics 
in the Toomre diagram. The metallicity gradient for 
thick disk stars with -0.9 < [Fe/H] < -0.3 is signifi- 
cant with [Fe/H] = -0.255|Z| -I- 0.02. RHB stars from 
the thick disk with \Z\ < 5 kpc have a peak metallic- 
ity of [Fe/H] ~ —0.6 dex, a peak rotation velocity of 
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Figure 15. R vs. \Z\ diagrams for the Halo I (upper panel) and 
Halo II (lower panel) sub-populations. 



Viot ~ 170 km s^^ , and a vertical scale height of |Z| ~ 1.2 
kpc. (2) The division between the thick disk and the halo 
seems to be quite clear in the [Fe/H] versus |Z| diagram 
as clearly shown in Figure 8, where the edge of the thick 
disk could be as high as \Z\ ^ 8 kpc at [Fe/H] ~ —0.6 
and reduces to \Z\ 2 kpc at [Fe/H] ^ —1.5. Mean- 
while, there is a detectable rotational gradient in the 
thick disk and a nearly constant Galactic rotation in the 
halo in the vertical direction to the Galactic plane. (3) 
Two halo sub-populations have been identified, the Halo 
I and the Halo 11, for stars peaking at [Fe/H] —1.3, 
based on the V^ot and the [Fe/H] — \Z\ diagram. The 
Halo I mainly clumping at i? < 10 kpc with a sign of 
metallicity gradient, while the Halo II clumping at two 
regions, both i? < 10 kpc and i? > 10 kpc, shows a 
negligible metallicity gradient. 

The above results are somewhat limited by the pres- 
ence of substantial scatter around these relations, which 
probably derives mainly from the large errors of the 
proper motions, and the limited sample of RHB stars 
with precise and accurate proper motions. Further study 
on a large sample of RHB stars with accurate proper mo- 
tions is desirable to clarify the two halo sub-populations. 
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In particular, the distance and proper-motion data for 
distant stars may be provided by the GAIA project. 
Combining these with the upcoming LAMOST spectro- 
scopic survey, we hope to identify more RHB stars and 
to analyze the RHB sample for structure in the metallic- 
ities, spatial coordinates, and Galactic velocities. 
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